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broken by the presence of the gluino mass. However, the renormalization of the supercurrent can be 
realized in a scheme that restores the continuum supersymmetric WTi (once the on-shell condition 
is imposed). The general procedure used to calculate the renormalization constants and mixing 
coefficients for the local supercurrent is presented. The supercurrent not only mixes with the gauge 
invariant operator T M . An extra mixing with other operators coming from the WTi appears. This 
extra mixing survives in the continuum limit in the off-shell regime and cancels out when the on-shell 
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I. INTRODUCTION 



SUSY or fermion-boson symmetry is one of the most 
exciting topics in field theory. From a theoretical point 
of view SUSY plays a fundamental role in string the- 
ory. There are many strong phenomenological motiva- 
tions for believing that SUSY is realized in Nature in a 
spontaneoulsy broken form. The SUSY breaking mecha- 
nisms are requested in order to produce a low energy 4- 
dimensional effective action with a residual N — 1 SUSY. 
For the other hand, non-perturbative studies of super- 
symmetric gauge theories turn out to have remarkably 
rich properties which are of great physical interest, as has 
been pointed out in {l}. For this reason, much effort has 
been dedicated to formulating a lattice version of super- 
symmetric theories (for a recent review in SUSY on the 
lattice with a complete list of relevant references, see Q). 
More recently, related interesting results in SUSY can be 
found inHflQEflflflyQy. Some of these for- 



mulations try to realize chiral gauge theories on the lat- 
ti„ with an exact chiral gauge symm e tty QQEQ. 
The lattice formalism is a powerful tool to extract non- 
perturbative dynamics of field theories and may be able 
to provide additional information and confirm or improve 
theoretical expectations. 

To formulate SUSY on the lattice we follow the ideas 
of Curci and Veneziano |17j . They propose to give up 
manifest SUSY on the lattice, and instead, to restore it 
in the continuum limit. In |l7j |. the Wilson formulation 
for the N = 1 SYM theory, which is the simplest SUSY 
gauge theory and corresponds to the SUSY gluodynam- 
ics, is adopted. For SU(iV c ) it has (N^ — l) gluons and the 
same number of massless Majorana fermions (gluinos), in 
the adjoint representation of the color group. 

SUSY is broken explicitly by the Wilson term and 
the finite lattice spacing. In addition, a soft breaking 
due to the introduction of the gluino mass is present. 
In [ljj it is proposed that SUSY can be recovered in 
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the continuum limit by tuning the bare gauge cou- 
pling, an d the gluino mass, m s - , to the SUSY point, 
rrinr. = 0, which also coincides with the chiral point. In 
EH HE 00, the DESY-Miinster-Roma Col- 
laboration have investigated these issues for the SU(2) 
gauge grou p (some first results have been obtained for 
ST/ (3) [25j), simulating the theory with a dynamical 
gluino using the multi-bosonic algorithm with a two- 
step variant called the (TSMB) algorithm |27j (while 
quenched results are in |28|V 

Another independent way to study the SUSY (chiral) 
limit in the Wilson formulation of Curci and Veneziano, 
is through the study of the SUSY WTi. On the lattice, 
it contain explicit SUSY breaking terms and the SUSY 
limit is defined to be the point in the parameter space 
where these breaking terms vanish and the SUSY WTi 
recovers its continuum form. These issues have been in- 
vestigated numerically in 0] in the on-shell regime. 

In this paper, the general procedure used to determine 
the renormalization constants and mixing coefficients for 
the local definition of the supercurrent, in the off-shell 
regime, is explained. It is shown that, when the opera- 
tor insertion involves elementary fields, the supercurrent 
not only mixes with the gauge invariant operator T^, as 
have been claimed in |17| . The supercurrent contains 
also non-Lorentz covariant terms which survive in the 
continuum, in the off-shell regime. These non-Lorentz 
breaking terms cancel out when the on-shell condition 
on the gluino is imposed and the continuum SUSY WTi 
is recovered. Preliminary studies have been presented in 
HQ. 

The paper is organized as follows. In Sec.|n]the Curci 
and Veneziano lattice formulation of the N = 1 SYM 
theory is presented, together with the lattice action and 
the vertices used for the calculation. In Scc. lIIII thc SUSY 
WTi on the lattice are written and the renormalization 
procedure explained. The calculation of the renormaliza- 
tion constant for the supercurrent is presented in Sec. lIVI 
Discussions and outlook are summarized in Sec. E] In 
Appendices ^ [5] an d 10 some details of the calculation 
are showed. 



II. LATTICE FORMULATION 

In the Wilson formulation of the N = 1 SYM theory 
0|, the gluonic part of the action is the standard pla- 
quette one 

S 9 = f EE f 1 - ^ReTriV(z)) , (2.1) 

where the plaquette operator is defined as [3^ 

P^{x) = Ut{x)Ul{x + P)U v {x + fiU^x) , (2.2) 

and the bare coupling is given by (5 = 2N c /g$. For Wil- 
son fermions the fermionic part of the action reads 

S /= 5> 4 Tr 



I 

2a 



A(z)(7m ~ r)UUx)X(x + afi)U^(x) 



X(x + a£)(7 M + r)U„{x)X{x)Ul{x) 
4r N 



rn 



A(a:)A(a!) 



(2.3) 



where uiq is the gluino bare mass and a is the lattice 
spacing. The fermionic field (gluino), X(x) — X a (x)T a , 
is a Majorana spinor in the adjoint representation of the 
gauge group. The symbol Tr implies the trace over the 
color indices. The normalization is given by Tr(T a T b ) = 
TjSai,. In this paper, only the case N c — 2 is considered, 
for which the adjoint gluino field is expressed in terms of 
Pauli matrices as 

(2.4) 

fc=i * 

The gluino field X(x) satisfies the Majorana condition 

X(x) = CX T (x), (2.5) 

where C = 7270 , is the charge conjugation operator. Our 
matrix convention for the Euclidean 7 matrices is as fol- 
low, 



A = E^ fc A* 



7o 




and 



7ft 




(2.6) 



(2.7) 
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The matrix 75 is defined to be 



75 = 71727370 



and the matrix a pv is 



1 

-1 



The anti-commutator property is 

{ipHv) = 26 ^ . 



(2.8) 



(2.9) 



(2.10) 



Finally, the Wilson parameter is fixed to be r = 1. 

SUSY is not realized on the lattice because, as the 
Poincare algebra, a sector of the superalgebra, is lost. 
SUSY is explicitly broken in the action (|2. 112.3(1 by the 
lattice itself, by the gluino mass term and by the Wilson 
term. Nevertheless, one can still define some transforma- 
tions that reduce to the continuum su per symmetric ones, 
in the limit a — > 0. One choice is [32I \?>?\ [5^ : 

6U ll (x) = -ag Q U p {x)ax)lpK x ) - 

ago£(x)jnX(x + ap,)U M (x) , 

SUl(x) = ag^{x) 1 ^\(x)Ul(x) + 

agoUl(x)£(x)^^\(x + afi) , 



5\{x) 



9o 



-<J P tG P t{x)£,{x) . 



SX(x) = —£,{x)a pT Qp T (x), 
9o 



(2.11) 



where £(x) and £(x) are infinitesimal Majorana fermionic 
parameters, while Q pT {x) is the clover plaquette operator, 

Q P r{x) = -^{P p r(x)-P TP {X) + 

P-p.-r(x) - P-r,-p{x) + Pr.-p(x) - 



P_p, T (x) + P-r,p{x) - P p ,-r(x)j . (2.12) 

Weak coupling perturbation theory is developed by writ- 
ing the link variable as 



U^x) = e- aA - 



aA lx (x) 



(2.13) 



and expanding it in terms of go. Here the gluon field is 
defined to be A M (x) = -ig Al\{x)T b . 

In order to calculate the 1-loop corrections to the 
SUSY WTi (which correspond to 0(#q)), we need two 



kinds of gluon-gluino interaction vertices. The gluon- 
gluino vertex, 



Kp C (P,<l) = 90 f 



abc 



p p a q p a 



. . ,Ppa g^a, 
ir sin — - — — — , 
V 2 2 ' 



(2.14) 



two-gluons-one-gluino vertex, 

v 2 t cd (k,p) = \agi{r e r hd + f ade r bc ) 

. ,Vpa q p a p p a q»a ] 

ilpSm{-^- + -^-)~rcos{-^- + -^-) 6 pv (2.15) 

and the three-gluons vertex 

G%% x {kMM) = *gof ab 

k v a . ,k 2 u,a ki p a 
6vx cos(— )sin(-| 1-) + 

r <ki\a . k v a k 2u a 
^ACOs(-2-)sin(— — ) + 

c ,k 2p a . k ix a k x a 
^cos(^-)sin(- — ) 



(2.16) 



These vertices are similar to QCD and the only differ- 
ence is that the fermion is a Majorana fermion in the 
adjoint representation of the gauge group instead of the 
fundalmental one. 

III. SUSY WTI ON THE LATTICE 

The vacuum expectation value of an operator O is de- 
fined to be 



(O) = J dUdXOe-^" 



(3.1) 



where S total is the total action on the lattice. By apply- 
ing an infinitesimal local supersymmetric transformation, 
with a localized transformation parameter £(x), the lat- 
tice WTi is written as Q, 

SO 



(OVMx))-2m (Ox(x)) 



6m 



O 



6S, 



GF 



5=o 



(OX s (x)) 



O 



6S 



FP 



6^) 



C=o 



(3.2) 



where Sgf is the gauge fixing term, Sfp is the Faddeev- 



Popov term and j|^j|£ =0 represents the contact terms 



4 



(see Appendices [X] an d [El f° r definitions) . This WTi is 
also discussed in |34|. Xs(x) is the symmetry breaking 
term coming from the fact that the action is not fully 
invariant under (|2.11|l . Usually X$(x) is a complicated 
function of the link variables and the fermionic variables 
[3^ . and its specific form depends on the choice of the 
lattice supercurrent. 

Let us define the lattice local supercurrent as 

2i 



S p (x) 



90 



Tr <^ G pT {x)a pT -i^X(x) 



while is the symmetric lattice derivative 
V M /(a;) = 2~(/(a; + afi) - f(x - afi)) 
and x( x ) corresponds to the gluino mass term 



X{x) 



(Jo 



Tr \ g pT (x)a pT X{x) 



(3.3) 



(3.4) 



(3.5) 



In order to renormalize the lattice WTi the operator mix- 
ing has to be taken into account. The standard way to 
renormalize the supercurrent is to define a substracted 
Xs, whose expectation value is forced to vanish in the 
limit a — > j^, [3] . In the case in which the operator 
insertion O in Eq. I|3.2[l is gauge invariant, Xs mixes with 
the following operators of equal or lower dimension [3] 

X s (x) = X s (x)-(Z s -l)VpS p (x) 

-2m X {x) - Z T V Mx) , (3.6) 



where the current T p reads 



Tp{x) = — - Tr \ Q^ v {x)^) v X{x) 



(3.7) 



On the other hand, if the operator insertion O is non- 
gauge invariant (i.e., the one involving elementary fiels), 
the gauge dependence implies that operator mixing with 
non-gauge invariant terms has to be taken into account 
in the renormalization procedure. In this case Eq. (|3.6() 
is modified as QE3 

X s (x) = X s (x)-(Z s -l)V fl S fi (x) 

-2rh X (x) - Z T X7 p T p {x) - Z B] B 3 . (3.8) 

i 

The Bj's denote the occurence of mixing, not only with 
non-gauge invariant operators but also mixing with gauge 



invariant operators which do not vanish in the off-shell 
regime (but vanish in the on-shell regime) . Consider, for 
example, the gauge invariant operator 

(3.9) 



B a = -Tryr P (D T g pT (x))X(x) 

which is zero imposing the equations of motion (thus, 
is not considered in |23j), but in the off-shell regime is 
non-zero and must be considered |38j. Other non-gauge 
invariant operators, which should be included in Bj are 

Si = -d p A p j9X,B 2 = -A p d p jdX, B 3 = - /Ad p d p X , 
9 9 9 

(3.10) 



(also reported in |32j). Finally, non-Lorentz covariant 
terms coming from V p S pi the gauge fixing term and con- 
tact terms, which appear in the off-shell regime, should 
also be taken in consideration. Because the Bj do not 



appear in the tree-level WTi, Z B] should be 0(g 2 ) 

Substituting (|3~%j) in (|3~2|) we obtain the renormalized 
WTi 

Z s (OVpSp{x)) + Z T (OVpT p {x)) - 
2(m - m)Z-\Ox R {x)) + Z CT ^i^| ?=0 

ZG F (o^k =0 )-z F Jo s -£^k= 



3 



sax 

Z Bj {OB j )=Q 



(3.11) 



The contact terms, Faddeev-Popov term and gauge fixing 
term should be renormalized, that is why in Eq. (|3.11|l 
the renormalization constants Zct, Zqf and Zpp are in- 
troduced. (OBj^j can in principle do mixing with S p and 
T p |29j. This implies that Sp not only mixes with T p as 
was predicted in [l7^. but extra mixing with gauge vari- 
ant operators and/or gauge invariant operators, which 
do not vanish in the off-shell regime, can appear. These 
extra mixing vanish by setting the renormalized gluino 
mass to zero and by imposing the on-shell condition on 
the gluino. 

In the continuum, the existence of a renormalized 
SUSY WTi 

d fl S* = 2m R Z xX , (3.12) 

is generally assumed, where Sr is the renormalized su- 
percurrent and mn is the renormalized gluino mass. For 
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to/?, = 0, we have SUSY while a non-vanishing value of 
?77a breaks SUSY softly. It is generally assumed that 
SUSY is not anomalous (Eq. 13.12JI holds) and only the 
mass term is responsible for a soft breaking. However, 
in |sd | the question of whether non-perturbative effects 
may cause a SUSY anomaly has been raised. 

It is tempting to associate the normalized continuum 
supercurrent as 



(3.13) 



in analogy with the lattice chiral WTi in QCD. This anal- 
ogy fails, as has been pointed out in |3t3 |. Explicit 1-loop 
calculation may shed some light on this issue. If the 
correctly normalized supercurrent coincides with l|3.13[) . 
then it is conserved when 777 ft = 0. This is the restoration 
of SUSY in the continuum limit j^ . 

By using general renormalization group arguments (see 
for example, one can show that Z$, Zt and Z x , 

being power substraction coefficients, do not depend on 
the renormalization scale /j, defining the renormalization 
operator in Eq. (|3.6|) . This imply that Zg = Zs(go, ra^a), 
Z T = Z T {g ,m a) and Z x = Z x (g ,m a). 

In this paper, we are interesting to calculate the renor- 
malization constant for the local supercurrent (|3.11l) and 
compare with Monte Carlo results in [2^. Notice that 
the relation between the 1-loop perturbative calculation 
and the numerical one is ZtZ^ 1 = Zt\i-Ioop- This is 
because, Z s = 1 + 0(g^), while Z T — 0(gl). So it 
is enough to calculate the coefficient Zt in 1-loop lat- 
tice perturbation theory (LPT). The numerical estimates 
are 2^] ZtZ^ 1 = —0.039(7) for the point-split current 
and ZtZq 1 — 0.185(7) for the local current, both at 
(3 = 2.3. An estimate of ZtZ^ 1 for the point-split cur- 
rent at f3 = 2.3 can be obtained from the 1-loop per- 
turbative calculation in [3^. At order g^ the value is 
Zt\i-Ioo P — —0.074 [3^]. In this paper, the calculation 
of Zt\i-Ioop for the local supercurrent is presented. 

In principle, each matrix element in Eq. H3.11J) is pro- 
portional to each element of the T-matrix base 

r = {1,75, 7a, 757a, <7 ap } , (3.14) 

but in order to determine Zt, it is enough to calculate in 



Eq. (|3.11(l the projections over two elements of the base 
(EH. 

IV. RENORMALIZATION CONSTANTS 

We are now considering each matrix element in Eq. I|3.11|) 
with O (a non-gauge invariant operator) given by 

0:=A b v (y)X a (z), (4.1) 

In Fourier transformation (FT) we choose p as the out- 
coming momemtun for the gluon field and q the in- 
coming momentum for the fermion field A (see Fig. ^| . 
Each matrix element can be written as 

(Al(y) \ a {z) C{x)) ^ D F (q) ■ (C(p, q)) amp -D B (p).S ab , 

(4.2) 

where (C(p,q)) amp can be, i.e., V^S^, V M T M , etc., with 
the external propagators amputated, Dp{q) and Db(j>) 
are the full fermion and gluon propagators, respectively, 
while S a b is the color structure, similar to all diagrams. 
The non-trivial part of the calculation is the determina- 
tion of (C(p, q))amp for each matrix element in Eq. H3.11JI . 
(Ox{x)): is not considered as we set the renormalized 
gluino mass to zero. 

In order to determine Zt one should pick up from each 
matrix element of Eq. (|3.11(l those terms which contains 
the same Lorentz structure as S p and T M , to tree- level. 
Those operators which do not contain the same tree-level 
Lorentz structure than S p and T p do not enter in the 
determination of Zt- Below, we present the tree-level 
values of the different operator of Eq. (|3.11() . The calcu- 
lation is straighforward. 

For the case of the supercurrent H3.3(l . the tree- level 
part reads 



d T A p (x))a pT ^^X(x) 



(4.3) 

Using Tr(T a T b ) = ^S ao for the traces and the antisym- 
metry of a pT this expression becomes 

SW{x) = -2S ab d p A b T (x)a pT ^X a (x) . (4.4) 

or in FT, 
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f ,4 f d*p f d*q 
= -2z5 ab a pTlfl JdxJ-^J-^ 

9-A f ^ f ^ v 

= -2zS ab a pT ^ J — 4 J ^ x 
xS(r-p+g)p P 4(p)A°(g), 



(4.5) 



7„ and |tr(7 Q 7 5 (V M S' A1 ) am p) the projection over 7^75 (tr 
is the trace over the gamma matrices which should not 
be confused with Tr, the trace over the color indices), it 
is easy to demonstrate that 

(0) 



1 



FT 



so we can define the vertex 



(4.6) 



Concerning the operator T M in Eq. (|3.7|l , the tree-level 



is 



H{p a Pv ~ Palv ~ P 2 5av + P ■ qS av ) (4.11) 

where jtr(7 M 7 p ) = S pp and 3^(757^7^7^) = e^, vpcr , 
while 

(o) 



jtr( 7^75 ( V,,S',, 



T^\x) = iS ab (d^A b T (x) - d T Al{x)) lT \ a (x) (4.7) 

after the FT, we define the corresponding tree-level ver- 
tex, 



Also 



^tr( 7a{ V M T M 



amp 



(0) 



amp 



FT 



2ip p q a e 



voLpa 



(4.12) 



FT 



T u,t(P> l) = S a b (jP$pr ~ Pplr) ■ 



(4.8) 



i{p a Pv - p a q v - P 2 8 av + P ■ q8 a „) (4.13) 



The tree-level expression for the amputated matrix el- 
ement (OV fiS^x)} (using the notation in Eq. 14.2fl is 



and 



(0) 



FT 



0. 



(4.14) 



V M 5 M 



(0) 



FT 



2(p- q) l xCr P ul l j,Pp, (4.9) 



Concerning the gauge fixing term, the tree-level value 
can be read from Eq. I|B9|) 



while the tree-level expression for the amputated matrix 
element (OV pT^x)) is 

(0) 



/ SS G f 1 



(0) 



FT 



ipp v 



(4.15) 



amp 



V T 



FT 



V7* 



' P ' qiv) ■ (4-10) and the projections are 



In our convention, = i(p — q)^ 1 is the momentum 
transfer of the operator insertion. 

From Eqs. (|4.9I4.10|I it is easy to see that, for p = q, 
a condition which would greatly simplify the calculation 
because implies that the operator insertion is at zero mo- 

/ N(0) / x(0) 

mentum, ( V A1 5 A1 (x) I = I V pT^x) I = 0. So the 

V / amp V / amp 

tree-level of and V^T^ can not be distinguished 

at zero momentum transfer. In order to determine Zt, 
different tree-level values of S 1 ^ and are needed. To 
differentiate these tree-level values, general external mo- 
menta, p and q are required. 

The value of the projections over 7 Q and 7 Q 7s for the 
different matrix elements in Eq. (13 . 1 1|> has been per- 
formed. Denoting |tr(7 Q (V A1 S l At ) amp ) the projection over 



1 / (5S GF 



(0) 



amp 



FT 



-2ip aVv (4.16) 



and 



1 / (SS GF , 



(0) 



amp 



FT 



0. 



(4.17) 



For the contact terms, the tree-level can be seen directly 
from Eq. (jB2j) (with a ->• 0) 



SO 
5£(x) 



(0) 



k=o > = 2i5{x-y) lu (\ a {y)\\z)) + 

5(x~y)(A:(y)a pT g b pT (z)) (4.18) 



or in FT (in the limit mo — > 0) 



2«7i/( — — L ]Sab - 2ip p —cr pi/ S ab . 

\i h) p 



(4.19) 
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The projections are 

1 / / SO . x(0) 

i tr V a {sW) k=a 



FT 



2i(p a qu - p ■ qS va + p 2 S a u) 



(4.20) 



and 



SO 



(0) 



amp 



FT 



-2i-p p q a e vap a ■ (4.21) 



Finally, the tree-level vertex for the operator in Eq. (|3.9() 
is 



{Bo)r b (P, q) = iSab(PpPr ~ P* S pT )l p 

while the projection is 



(4.22) 



1 



tr la (B ) 



(0) \ FT . 



amp 



i(j>uPv - P S av ) . (4.23) 



For the operators in Eq. i|3.10[l we have 

I FT 

amp) ^ iPisQa ? 

j FT 

amp) ^ IQvQa •> 

-tr( 7Q (B 3 ) amp) O a u 



and 



(4.24) 



itr( 7a75 (5 , 1)2 ,3) H)^0. (4.25) 

The renormalization constants can be written as a 
power of go 



Z 



_ 7 (o) 

operator — ^operator 



2 7 (2) 



operator 



(4.26) 



and also for the operators, a similar expansion can be 
done it 



(Operator} = (Operator)^ + g$ {Operator) ^ + ■ • • , 

(4.27) 

(2) 

where (Operator) , is the 1-loop correction while 
(Operator) , is the tree-level value. 

Substituting Eq. and Eq. !j4"27|) into Eq. l(3~TT|) . 

to order g"^, we obtain 



(l + . 9o 2 4 2) )(<OV M 5 M (x)) (0) + 

. 9o 2 (0V^(*)) (2) 



5o 2 4 2) (OV M T M ( a 



AO) 



so 



ss, 



(0) 



(2) 



9o (O 



GF , 



l€=o 



(2) 



(0) 



2 y(2) / n ^FP , 



(4.28) 



At tree-level we have, Z ( ° ] = 1,Z^ = 0,Z ( ° ] T = 
1, Z^l = 1, Zf } p = 0, = 0, so the lattice WTi is 



(0) 



(0) 



0. 



(4.29) 



which holds in our lattice calculation. Eq. Q4.29|l was 
previously determined in the continuum 40] . To order 
<7q the lattice WTi is 

{OV,S,(x)) {2) + 4 2) (OV M ^(x)) (0) + 

zP(ovMx)) (0) + 



so 



(2) 



z 



(2) 
CT 

(0) 



/ SO 



(0) 



(0) 



(2) 



<5Sj; 



(4.30) 



(0) 



Notice that the Faddeev-Popov term | g^o ) m 

Eq. (|4.28() . is already 0(g 2 ) (see Appendix lB|l and does 

not contribute to 1-loop order. In Fig. ^ the Feynman 

/ \ (2) 

diagrams for (OV,S m (j:)) (2) and /o^| f=0 \ are 

shown, while in Fig. [S] the non-zero contribution to con- 
tact terms are presented. 

Let us substitute the tree-level values of the operators 
in Eq. (|4.30|l using the projections over 7cn 

(2) 



amp 



^tr( 7a { 
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Zf^iipaPv - p a q v - p 2 8 av + p ■ q5 av ) + 
Z^'i{p a p u - p a q v - p 2 5 au + p ■ qS au ) + 



SO 



\£=o 



(2) 



amp 



Zfc2i(p a q v - p ■ q5 av + p 2 5 aiJ ) 



rtr 7. 



SS, 



GF , 



1 



-Z^( la OB 3 r=Q 



k=o 



(0) 



(2) 



(4.31) 



and the projections over 7 Q 7s, 
1 



(2) 



amp 



^tr ^75 [V^S, 

Z^-Jp c lippq a £ ua p a -\- 



4 tr l 7a75 l^) l? =° 

1. / / $Sqp 



(2) 

amp 
(2) 



+ 



j^tr<7 Q 7 5 O^-> w =0. 



(0) 



(4.32) 



(2") 

Our claim is that, in order to calculate Z\, we 



can substitute i^tr^O^) 



AO) 



p 2 5 av ) + Z^[ip v q a + Z^iq v q a + Z^iq 2 5 av 



(2) 



'(2)- 2; 



and 



i^ 2 i ) Tr(7 Q 7 5 OB i ) 1 "' ->■ 0, where Z Bj correspond to 
the renormalization constant in Eq. (|4. 2314. 2414. 25jl . No 
other are needed, because there are no other Bj's 
that would contribute with the same Lorentz structures 
appearing in the tree-level of Eqs. 14.31l4.32fl . 

Each matrix element in Eqs. (|4. 3114. 3211 has been cal- 
culated for general p and q (off-shell regime). To deal 
with the IR divergencies and rinormalize to 1-loop order 
the Kawai procedure is used , with the help of tabu- 
lated results in 0,^3 . O nce p and q have been extracted 
from the propagators through the Kawai procedure, the 
rest of the integral depend on the loop momenta which 
is numerically integrated. A similar renormalization pro- 
cedure has been used to calculate the 3-loop beta func- 
tion in QCD with Wilson fermions (44| and the 3-loop 
free-energy in QCD with Wilson fermions (for 
a complete study of the off-shell WTi in QCD see |47|). 



A°) 



Tipically, each matrix element contains f» 1000 terms (in 
particular dilogarithms functions depending on both ex- 
ternal momenta which coming from the diagrams with 
three propagators in Fig. QJ. After the numerical inte- 
gration, one can simplify the results in order to read the 
value of Zt by setting 



q and p ■ q = , 



(4.33) 



(see Appendix [UJ. This is still an off-shell condition (be- 
cause even if p 2 — q 2 , there are no other condition on 
this expression, i.e., q 2 = 0), but drastically reduces the 
number and difficulty of the expressions (for example, the 
dilogarithm terms simplify). 

Let us introduce, for simplicity, the notation, A = 
0V^(:r) + ^|£=o-0f^| e =o. Using Eq. we 

(2) 

get the following dependence on p and q for tr^7 a A) 

(2) 

and tr(7a7 5 A) , 

tr( 7a A) (2) m A 

lQ 2 PaPv + A 2 q 2 p a q v + 
(A 3 + M 3 )q 2 6 al/ + M X q 2 p v q a + M 2 q 2 q a q u + 
Piq 2 plS, a + P 2 q 2 qX a + ••• (4.34) 



and 



tr(7 Q 7 5 A) 



(2) FT 



A i q 2 p p q IJ e v 



(4.35) 



where the dots in Eq. H4.34J1 indicate that because the 
scmplification in Eq. I|4.33|) is used, some momenta de- 
pendence are missing or mixed with others, i.e., p ■ q5 va 
does not appear, while p 2 5 ua , is mixed with q 2 8 va , (see 
Appendix El for notation) . 

It is also interesting to see the Lorentz structure of the 
supercurrent, 

tr(-f a S^) (2) ^> N 1 qp tl p u p a + N 2 qq tl p v p a + 
Nzqpp,q v p a + N^qq^qupa + N^qp^qa + 
Neqq^qa + N 7 qppq u q a + N 8 qqpq v q a + 

QiqpaSpv + Q 2 q q a Su» + Qzq 

Qiqq v 5p a + Q^qP^ua + Q 6 q 

RiqPfiS^a + R.2qq^nva H (4.36) 

where the coefficient Ai, Mj, Qk, are tipically of the form 

(C n + C m Ln(aV)) (4.37) 
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while Pi,Nj,Rk, do not contain Ln(a q ) terms. Here, 
C n are lattice constants or numbers coming from the nu- 
merical integration and C m are rational numbers com- 
ing from the Kawai procedure. Notice that the Lorentz 
structures multiplying P t , in Eq. 1(4. 3414. 36|) are non- 
Lorentz covariant, even in the continuum limit (a — > 0). 

From Eqs. I(4.31I4.32I4.34I4.35(> the following condi- 
tions can be derived 



Ai = -2iZ^ ] - iZ { S> 



(2) 



A 3 + M 3 = 2iZ 



(2) 

s 



iZ^ ] - iZ 



(2) 



Mi 



-iZ 



(2) 
Si ' 



M 2 =-iZ B 



(2) 



(4.38) 



and 



An 



2iZ 



(2) 



iZ, 



(2) 



2iZ, 



2) 
CT > 



A 4 = -2iZ 



(2) 



2iZ. 



CT ■ 



(4.39) 

(2) 

The last two conditions can be explicitly solved for Z T , 



Z. 



(2) 



-iA 2 — iAi 



(4.40) 



Eq. (|4.4()(l is the only possible solution of the system 



(2) I 

7" 1 |1 — loop 



= 0.664. 



(I4.3dl4.3^ll for Z ( T Z> . Our result is Z, 
A VEGAS Monte Carlo routine to perform the 1 loop 
integration with 200 millions of points, using the GNU 
Scientific Library is used. To estimate the error we 
take the value given by the program which is ss 10 -5 for 
each integral. The calculation, once p and q has been ex- 
tracted from the propagators, involves around 1300 dif- 
ferent 1-loop integrals. For each diagram tipically we 
have 100 different integrals. That means that the error 
is around 10~ 3 . 

Let us compare our perturbative result with the nu- 
merical one 0, Z£ UM = Z T /Z S = 0.185(7). One 
has to observe that the definition used here for S u is 



not the same as in 
tNUM _ 1 7PT ' 



It is easy to demonstrate that 
Z^uiu _ Lz^ 1 [5fJ. To compare with the numerical re- 
sults one as to divide the perturbative value by 2 which 



gives, 



7PT 



1 7(2)1 

2 Z T ll- 



loop 



0.332. We are currently in- 



creasing the precision of the numerical integration to 400 
milion of points. A detailed presentation of the results in 
Eqs. I|4. 3414. 351) together with the result of each diagram 
is under way pi). 



V. DISCUSSION AND CONCLUSIONS 

In this paper, the SUSY WTi in 1-loop LPT has been 
investigated. A general procedure in order to get the 
renormalization constant for the supercurrent has been 
presented. In LPT it is possible to determine the value 
of the renormalization constant for the supercurrent from 
the off-shell regime of the SUSY WTi. The computation 
of each matrix elements of the WTi has been carried 
out using the symbolic language Mathematica. The pro- 
grams were completely wroted by the author together 
with the numerical code used for the integration. All 
the contributions have been calculated in the off-shell 
regime, and in order to get the value of the renormaliza- 
tion constant, a simplification in the external momenta 
(which still keeps the off-shell regime) has been applyed. 
We are currently increasing the precision of the numeri- 
cal integration and a detailed presentation of the results 
is the subject of a forthcoming paper [5^]. A reason- 
able good agreement of our perturbative result for the 
renormalization constant, ZJp T = 0.332, in comparison 
with the numerical one, Z^ UM = 0.185(7), has been 
achieved, taking in consideration the fact that in the nu- 
merical simulation = 4/2.3, which still corresponds to 
the non-perturbative region. We observe that, at least 
at 1-loop order in perturbation theory, Zt is finite. This 
result may have some theoretical implications which we 
are currently investigating. Also, the determination of 
Zs, using another kind of gamma projections is under 
investigation. It would be interesting to calculate Zs in 
order to check the trace anomaly and the exact renormal- 
ization expression for Eq. 1)3.13(1 . An important point to 
stress here is that, even in the continuum limit, we ob- 
serve in Eq. 14.34J1 Lorentz breaking terms which coming 
from the fact that we substituted X$ by the Eq. (|3.8|l . 
It would be interesting to see whether Eq. H4.34|l is the 
continuum off-shell WTi. The nice point is that, once the 
Zt has been determined, we can impose the on-shell con- 
dition on the gluino mass. The Lorentz breaking terms 
cancel out from Eq. 1)4.34(1 and the continuum WTi is 
recovered. At least to 1-loop order, we do not observe a 
SUSY anomaly in TV = 1 SYM, altough a more carefully 
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study is required. 
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APPENDIX A: PERTURB ATI VE CALCULATION 



In this appendix we follow the lines of |29U4(l ]. The lattice 
SUSY transformations of the gauge field (x) are not 
equal to the continuum ones. On the lattice the trans- 
formation of the gauge link U^ix) determines the trans- 
formation properties of A^(x). Writing the link variable 
as 



17 (a) = e - aA f ( x +##) , 



(Al) 



for the SUSY transformations of the gauge link we use 
the symmetric choice 

-ag t;(x + afi)^f^\{x + afyU^x) . 

These two equations determine the transformation be- 
havior of the field A^x) 29]. The FT for the gauge field 
is defined as the usual way 



A b (x) = / d A kA b Jk)e ik < x+ % ft ) 



(A2) 



Collecting all terms until order g$ we can write down the 
variation of the gauge field A b (x) as 

5Al(x) = i (^{x)^X b (x) + £(x + afi)^X b (x + afi^j + 

-^agofabc ( £(x)7 M A c (x) - f (x + a/i)7 M A c (x + afi) ) A a - 



i 

24 



a 2 gl ( 25 a bScd ~ S ac Sbd — 8 ad 5bc )A c li A lt x 



x i (xhp.X a (x) + £(x + a/i)7 M A a (x + a/t) 



(A3) 



which reduces to the continuum SUSY transformation 
5A"(x) = 2i^7 AI A a (x) in the continuum limit a — * 0. 
Because in this paper we fix N c — 2, some semplifications 
appear 



1 



Tr ^ T a T „ T c T ay = -(S ab 6 cd - f abe fcde), 

fabefcde = {S aC S b d ~ S ad S bc ) . (A4) 



Using the Eq. (|A3|) , it is possible to determine the dif- 
ferents pieces of the WTi in Eq. H3.11JI . i-e., the con- 
tact terms, the gauge fixing term and the Faddeev-Popov 
term. They are necessary in order to calculate the Feyn- 
man rules for 1-loop order calculation. In appendix [5] 
the vertices coming from these pieces are presented to- 
gether with the ones coming from the supercurrent. 

APPENDIX B: VERTICES 

Let us determine the contact terms, lg=0- First 
of all, the variation of the operator insertion, O = 
A a v (y)X b (z), is 



50 = 5Al(y)\ b {z) + Al{y)8\ b {z) 



(Bl) 



Substituting I|A3J| into l|Bljl . after some algebra, we ob- 
tain 

SO , 



\S=a 



i5(x-y)~fJ\ a (y)\ b (z)) + 



iS(x -y- az>)7 y ^A a (2/ + av)\ b {z)j + 
ag fdacS(x - y) ly (\ c {y)Ai(y)\ b (z) ^ 
ag fda C 5(x -y- az>)7„( X c (y + av)A d v (y)\ b (z) 



24 



a g Q 25 ae 5 cd - 5 ec 5 ad - 5 ed 5 ac 5(x - y) x 



x (AUy)Ai(y)X e (y)X t '(z: 



--aVo ( 26 ae S cd - S ec 6 ad - S ed 6 ac )S(x-y- av) x 
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x (A c v (y)AUy)\%y + aP)\<>(z) + 



5{x-y)(Al{y)a pT Q b pT {z)y (B2) 
where £ pT (z) - -ig g b pr {z)T b . 

The part of the lattice action corresponding to the gauge 
fixing is defined as 



X p 

-^4 



(B3) 



= y E (E v p qc ^pW) (E V ' QC ^)) 



a: P 



where V b p ack f(x) = i f /(x) — f[x — ap)\ . The variation 
of the gauge fixing term (|B3f) can be written as 



X p 

(B4) 

This results in the contribution of the gauge fixing term 
into the WTi as 



^7 P (\ c {x)V s p orw V b T ack iyA c T (x) + 

A c T {x-ap)^A a M\ b {z) 
l -ag Q f ecllp (\f{x) ^A;wf p orw W b T ack A c T (x) - 
A e p (x - ap)Vf,° rw V b r k AUx - ap)]Al{y)X b {z) 



.9o7p^2(5 ec (5/d - SefScd - SedSfcj x 

x (^X e (x) (^AlA d p V f p orw V b T ack A c T {x) + 
Af(x - ap)A%x - ap)Vf°™V b T ack A c T (x - ap)) x 



xAl(y)\ b (z] 



(B5) 



where V{,° rw f(x) = I (/(a: - a/5) - f(x)). 

Finally, the expansion of the Faddeev Popov action can 

be written as 



« 2 E 



x,p<0 



fi a {x)[5 ab + -ag Q Air b + 



1 

12 



^ a (x)((5 Qb --a ffo yl p r cfc + 



12 



(B6) 



and the contribution of the Faddeev-Popov term into the 
WTi is 



-^/-E7p 



^(a;)A c (x) (v h (x) + f] h (x + apU + 
ff(x - ap)X c (x) [r) h [x - ap) + r) h (x)J - 
fj 9 (x + ap)X c (x) \r) h {x + ap) + r) h (x)j + 



fj g (x)\ c (x)(ri h (x)+ri h (x~ap) 



A a Ay)\\z) 



l 90 jgch jdcf ^ j p 



fj s (x)X f (x)A d (^ h {x) + r) h {x + ap)^j - 
ff (x - ap)X f (x)A d (x - ap) \r) h (x - ap) + r) h (x) 
rj 9 (x + ap)X f {x)A d L h (x + ap) + V h (x)^ + 
fj 9 (x)X f (x)A d (x - ap) (r} h (x) + r) h (x- ap) 



xA a „(y)\ b (z] 



12 



fj°(x)(\ c (x)A d p + A p X d (x)^ 

r} h (x) - T) h (x + ap)j + 
fj 9 (x-ap)(x c {x)A d p (x-ap) + 

A c p (x - ap)X d (x)^j [r, h {x - ap) - V h (x)^j + 
f\ a {x + ap) (x c {x)A%x) + A p (x)X d (x) 
r] h {x + ap)-T] h {x)) + 



12 



7] 9 {x) X c (x)A d p (x - ap) + A c p (x - ap)X d (x) 



r] h (x) — rj h (x — ap) 



A a M\»{z) , (B7) 



where A c p = A c p (x + °^-). It is possible to calculate the 
vertices and the corresponding Feynman diagrams, up to 
order g$, from Eq. jHUBUBIIl in FT. 

Regarding the contact terms in Eq. I|B2|) . all the con- 
tributions to order <7q are zero except for the last line 
of Eq. I|B2|) . The corresponding non-zero Feynman dia- 
grams are shown in Fig [5] The vertices used here are, 
the two-gluons vertex, 



1 



cos( 



P T 



sin( 
sin( 



ki p a 




2 


2 


ki T a 


k 2rO. 


2 


2 


ki p a 


k 2p a 


2 


2 


k\ T a 


k 2T a 


2 H 


2 



' 2 ' 

,kipa 
2 



)sin(^)cos(^) 



, k 2p a 

k\ T a^ 
2~ 



)sin(-y-)cos(-|-) 



VL^rSml^ t~ ~z ) [ sm(fc lct a) 



sin(fc 2a a) 



(B8) 



and the three-gluons vertex, which we do not reported 
here and gives a zero contrbution to the last diagram of 
Fig.H 

For the gauge fixing terms in Eq. |B5j), we need the 
vertex with one-gluon-one-gluino (which is similar to 
Eq. 14.15f) in the continuum limit, 

gf i P t(p> k ) = --^Sablp sin(fcpa) sin(^) , (B9) 

the vertex with two-gluons-one-gluino 

— — f ce 

GF 2Tp (p 7 q,k) = 

2ffo , / . ,k p a q p a . q p a . q T a 
— /ec/|7psm(— + — )sm(— )sm(— ) - 

7 r sm 1 ) sin ) sin — — ) 

it y 2 2 ' y 2 ' V 2 ' 



(BIO) 



and finally the three-gluons-one-gluino vertex (non- 
symmetrized) 

efdc 

GF 3p<7T {p,k,q,t) = 



— gflo \25ec5fd - SefScd ~ 3ed$fcj Jf>Spa X 

xsm(-|- + ^ + ^-)sin(^-)sm(^-). (Bll) 



For the Faddeev-Popov terms in Eq. (|B7|I . we need 
one-gluino-ghost-antighost vertex 

cgh , 

FP p (p,-q,k) = 

4 50 rnchSr^ ( k P a \ ■ l Q P a \ l k P a 1p a \fT3TO\ 

—f 1^ 1p cos (-§-) sm(y) cos(-^ p —) (B12) 

p 

and one-gluino-one-gluon-ghost-antighost vertex 

- — ^cdqh 

FP lp (p,t,-q,k) = 



2i 



gl{f 9ce f edh + f 9de f ech )i P 



sm(y) sm(y) cos (^~ + ~y 



q P a . 

2 1 



(B13) 



As we can see from Eq. <|B12|) and <|B13ll the vertices 
are already order g and g^, so pluging into Eq (|4.28|) is 
already more than O(<7o)- This imply that the Faddeev- 
Popov terms do not contribute to order g$ . 

Concerning the vertices of 5 M for a 1-loop calcula- 
tion we need the vertices corresponding to one-gluon-one- 
gluino, the two-gluons-one-gluino and finally the three- 
gluons-one gluino. They can be calculated from ((3.3(1 . 
The vertex one-gluon-one-gluino (using Eq. I|4.5|) 'l is 

Si^pMiP) = -— Sab<y p r7p cos(^) sin(p p a) , (B14) 

which reduce to the continuum one in the limit a — > 
(see Eq. 1)4.6(1 1. while the vertex two-gluons-one-gluino is 



1 



Ti-x -CI Tin 

cos( 



cos(^— + — — ) cos(— — ) cos(^- +pi p a) 



sin( 
sin( 



1 2 


f 2 


,P\ T a 


P2tCI 


- 2 


2 


Pi P a 


P2 P a 


2 


2 


p lT a 


P2 T a 


2 


2 



) cos(-^-) cos(-^- + p2 T a) + 



)sin(-^-)cos(— -) 



x . ,P2 T a ,P2 P a, 
) sm (^-) cos (^) 



j p t 



Out1 p sin( 



P\tCl p 2T a, 



) sin(pi a a) 



sin(p 2 ao) 



(B15) 



We do not presented here the three-gluons-one-gluino 
vertex because its contribution to the last Feynman di- 
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agram for the supercurrent, in Fig. ^ is zero by color 
considerations. 



others scalars, 



(C5) 



APPENDIX C: OFF-SHELL REGIME 

In order to separate the contribution of and at 
tree-level, we can not impose p — q, that would greatly 
simplify the calculation. We are forced to use general 
external momenta p and q (while the momentum trans- 
fer of the operator insertion is (p — q) 7^ 0, see Fig. [[). 
Once the external momenta has been extracted from the 
propagators, in order to get the value of Zt, the sem- 
plification p 2 — q 2 and p ■ q — is used. This is still an 
off-shell regime which simplify the dilogarithm functions. 

At 1-loop order, two propagators integrals are tabu- 
lated in |4l|, |42j while three propagator integrals on the 
lattice are tabulated in in terms of lattice constants 
plus the following continuum conterparts 



Io;l f i;2 f j, V ;3 f tu P (p,q) = — / d k 



1, A^, kpk Vl k^kykp 

'' k 2 (k+ P ) 2 {k + q y 



(Cl) 

With the help of |48|, |49j one can give the expression for 
Io[p,q) and write down recursively h^,{p,q), I 2f j,v(p,q), 
I^vpiPiq) in terms of the scalar functions p 2 , q 2 , p ■ q 



and Iq , plus Lorentz structures. As an example [4£ 
Io(p, q) is a complicated function of p and q, in terms of 
the dilogarithm as follows 



Io(p,q) = 



Li-, 



p-q-A 



— Lis 



p-q + A 



1 (p-q~A\ f{q-p) 



-Ln 



Ln 



2 "\p-q + AJ \ q 
where A is the triangle function defined as 



and 



A 



Li 2 (x) 



(p-qY 



2 2 

p q 



Lnt 
t-l~ 



dt 



(C2) 



(C3) 



(C4) 



is the dilogarithm. 

Following the reference 0] where a tensor decomposi- 
tion of hftip, q), hfiuip, q), h^upip, q) is used, it is shown 
that all the integrals can be written in terms of Iq and 



where 



h(p,q) = 



2 T 1 (y-p) 2 \ t A^-p) 2 

q Ln - - p ■ qLn 



q 2 p -(q-p) 



(C6) 



The integral is symmetric in [i and v as well as 
under p <-» q and hence has the following tensor decom- 
position 

S -fpAl B (p,q ) + 



hp, v = 5,0,1a + yP^Pu 

p^qv + q^Pv - 
5 



-yp -q)ic 



q^q v 



)lB(q,p) , 



(C7) 



where I a, Ib and Ic are symmetric under p «-> q and 
tabulated in [4^ . In this reference an explicit expression 
for i^i/A is presented, which is quite complicated and we 
do not reported here. 

The general result for arbitrary p and q using l|C2() . 
(IC6I) . (|C7(1 and the corresponding expression for 
(in |49() contains a huge quantities or terms (sometimes 
up to 1000 terms). Therefore a semplification which still 
leave us in the off-shell regime is required. Let us rewrite 
(|C3|> in the following way 

(p ■ q? 



A 2 



2 2 
-p q 



p 2 q 2 



1 



(C8) 



where p ■ q = pq cos a, where < a < n. This imply that 

(p — q) 2 — p 2 + q 2 — 2pq cos a - 

By using Eq. I|C8J) it is possible to simplify I , I\ J 2 
and I3. In fact, 



(C9) 




cos cr 



Substituting 

in (|C2(I we have 
1 



ip gVl — cos a 2 
p q cos a — ip 



Li 2 



qVT 



cos a z 
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Li 



. fpq cos a + ip qy/l — cos a 2 



1 fpq cos a — ip qyi — cos a 2 
-Ln = I x 

2 \pqcosa + ipqvl — cos a 2 , 



xLn 

simplifying we have 
1 



P 



q 2 — 2p q cos a 



(CIO) 



Jo 



ip q sin a 



Li 2 - (cos a — i sin a) 
9 



Li 2 ( - (cos a + i sin a) 



1 / cos a — i sin a\ f p 2p 
— Ln Ln — + 1 cos a 

2 V cos a + i sin a / V <jr q 



and finally 



(Cll) 



Jo 



ip q sin a 



Li 2 (^exp- M )-Li 2 (^exp"*; 



2 P 

1 cos a 



, (C12) 



or 



1 



p q sm a 
aLn( — 



i(Li 2 (^ex P — ) 
l V q 



Li 2 (^exp lQ ; 



1 



2p 



(C13) 



Using Eq. iCTl we can now simplify the recursive ex- 
pressions for Ii , 7 2 and .Z3 . Let us define 

P/i 



(C14) 



where clearly = 1. The semplification in Eq. 
correpond to a = % (then, we have cos a = and sin a = 
1) and p 2 = q 2 , which correspond to the substitution 
p — > q in all the results 
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FIG. 1: Diagrams contributing for the supercurrent and the 
gauge Sxing term. The grey blob correspond to the operator 
insertion in which Bows a momentum (p — q). 




